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ABSTRACT 



Context. The final performance of current and future instruments dedicated to exoplanet detection and characterization (such as 
SPHERE on the European Very Large Telescope, GPI on Gemini North, or future instruments on Extremely Large Telescopes) is 
limited by uncorrected quasi-static aberrations. These aberrations create long-lived speckles in the scientific image plane, which can 
easily be mistaken for planets. 

Aims. Common adaptive optics systems require dedicated components to perform wave-front analysis. The ultimate wave-front 
measurement performance is thus limited by the unavoidable differential aberrations between the wavefront sensor and the scientific 
camera. To reach the level of detectivity required by high-contrast imaging, these differential aberrations must be estimated and 
compensated for. In this paper, we characterize and experimentally validate a wave-front sensing method that relies on focal-plane 
data. 

Methods. Our method, called COFFEE (for COronagraphic Focal-plane wave-Front Estimation for Exoplanet detection), is based 
on a Bayesian approach, and it consists in an extension of phase diversity to high-contrast imaging. It estimates the differential 
aberrations using only two focal-plane coronagraphic images recorded from the scientific camera itself. 

Results. We first present a thorough characterization of COFFEE'S performance by means of numerical simulations. This charac- 
terization is then compared with an experimental validation of COFFEE using an in-house adaptive optics bench and an apodized 
Roddier & Roddier phase mask coronagraph. An excellent match between experimental results and the theoretical study is found. 
Lastly, we present a preliminary validation of COFFEE'S ability to compensate for the aberrations upstream of a coronagraph. 

Keywords, instrumentation: adaptive optics, instrumentation: high angular resolution, techniques: image processing, meth- 
ods: numerical, methods: laboratory, telescopes 



O 1 ■ Introduction et al. 2010) to compensate for these NCPA. 

■ Several techniques dedicated to high-contrast imaging system 

Exoplanet imaging is one of the main challenges in today's as- p tim i zation have been pr0 p 0se d for future systems. Some of 

^ tronomy. A direct observation of these planets can provide infor- them rely on a dedicated wave-front sensing hardware (Guyon 

^ mation on both the chemical composition of their atmospheres et al 2009X others use scientiflc foca i plane data assuming small 

. . and their temperatures. Such observations have recently been aberra tions. Speckle nulling iterative techniques (Borde & Traub 

m £ piade possible (Kalas et al. 2008; Marois et al. 2008; Lagrange 2 006; Give'on et al. 2007) estimate the electric field in the de- 

k> et al. 2009), but only thanks to their high mass or their wide ap- tector plane using at least three images The tec hnique proposed 

^ pare nt distance from their host star. by Baudo z et al. (2006) relies on a modification of the imaging 

& Being able to image an object as faint as an extra-solar planet system , but requires only one image. These techniques aim at 

very close to its parent star requires the use of extreme AO minimizing the energy in a chosen area ("Dark Hole"), leading 

(XAO) systems coupled to a high-contrast imaging technique, t0 a contrast optimization on the detector (Trauger et al. 2010; 

such as coronagraphy. Instruments dedicated to exoplanet imag- Baudoz et al. 2012) in a closed loop process, 

ing using these two techniques (SPHERE on the VLT, (Beuzit We have rece ntly proposed a focal-plane wave-front sensor, 

et al. 2007), GPI on Gemini North, (Macintosh et al. 2008)) are COFFEE (Sauvage et al. 2012), which is an extension of conven- 

currently being integrated. The performance of such systems tional phase diversity (Mugnier et al. 2006) to a coronagraphic 

is limited by residual speckles on the detector. These speck- system since COFFEE uses focal-plane images, it is possible 

les, which originate in quasi-static non common path aberrations t0 charact eri Z e the whole bench without any differential aberra- 

(NCPA), strongly decrease the extinction provided by the coro- tion Xnis method requires only two focal-plane images to es- 

nagraph and can be difficult to distinguish from an exoplanet. timate the aber rations upstream of the coronagraph without any 

To achieve the ultimate system performance, these aberrations modification of the coronagraphic imaging system or assuming 

must be measured and compensated for. The current-generation small aber rations. COFFEE'S principle and its application to the 

instruments, SPHERE and GPI, respectively rely on phase diver- apod ized Roddier & Roddier phase mask (ARPM) are described 
sity (Gonsalvez 1982) and an interferometry approach (Wallace 
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Fig. 1. Coronagraphic imaging instrument: principle 

in Section 2. In Section 3, we evaluate the quality of NCPA 
estimation by realistic simulations. In Section 4, we present 
the experimental results from the laboratory demonstration of 
COFFEE on an in-house adaptive optics bench (BOA) with an 
ARPM. Section 5 concludes the paper. 

2. COFFEE: principle 

2. 1. Extension of phase diversity to coronagraphic images 

Figure 1 describes the coronagraphic imaging scheme consid- 
ered in this paper. We consider four successive planes denoted 
A (circular entrance pupil of diameter D u ), B (coronagraphic 
focal plane), C (Lyot Stop), and D (detector plane). The optical 
aberrations are considered as static and introduced in the pupil 
planes A and C. The coronagraphic device is composed of a fo- 
cal plane mask located in plane B and a Lyot Stop situated in 
plane C. No particular assumption is made on the pupil shape or 
intensity. Thus, the description of COFFEE is compatible with 
several coronagraphic devices. COFFEE uses two images, i f c 
and i A c , recorded on the detector (plane D in Figure 1) that, as 
in phase diversity, differ from a known aberration, 4> div , to es- 
timate aberrations both upstream (4> u ) and downstream (4> d ) of 
the coronagraph. 

Considering the calibration of the instrument with an unre- 
solved object, we use the following imaging model: 



•foe 



a.h det *h c (cj) u ,cj) d ) +n 



foe 



i™ = a.h det * h c ((f> u + (f> div , 4> d ) + n 



div 



(i) 



where a is the incoming flux, h c the coronagraphic "point spread 
function" (PSF) of the instrument (i.e. the response of a corona- 
graphic imaging system to a point source), h det the known detec- 
tor PSF, n toc and n dlv are the measurement noises, j3 is a uniform 
background (offset), and * denotes the discrete convolution oper- 
ation. Such an imaging model can be used for any coronagraphic 
PSF expression h c . The measurement noises n foc and n dlv com- 
prise both photon and detector noises. Because calibration is 
assumed to be performed with high flux levels, we adopt a non- 
stationary white Gaussian model, which is a good approximation 
of a mix of photon and detector noises. Its variance is the sum of 
the photon and detector noise variances: er^[t] = <Xp h [i] + tr| et 
(Mugnier et al. 2004), with t the pixel position in the detector 
plane. The former can be estimated as the image itself thresh- 
olded to positive values, and the latter can be calibrated prior to 
the observations. 

We adopt a maximum a posteriori (MAP) approach and es- 
timate the aberrations <$> u and <f> d , the flux a, and the background 
/3 that minimize the neg-log-likelihood of the data, potentially 
penalized by regularization terms on <p u and <p d designed to en- 
force smoothness of the sought phases: 



(a,j3,4) u ,4> d ) 



argmin J(a, f3 7 (f> u , <f> d ) 



(2) 



where 



J{a,(3,(t> u ,4>d) = 2 



-{a.h A *h c (4> u ,4> d )+p) 



r foc 



if -(a.h d *h c (<t> u + <l> div) ct> d ) + l3) 



(3) 



where ||cc|| denotes the sum of squared pixel values of map x, 
cr f ° c , and cr A ™ are the noise standard deviation maps of each im- 
age, and 1Z is a regularization metric for the phase. 
Any aberration <fr 1S expanded on a basis {Zk} that is typically 
either Zernike polynomials or the pixel indicator functions in the 
corresponding pupil plane: <fi = J^k a kZk where the summa- 
tion is, in practice, limited to the number of coefficients consid- 
ered sufficient to correctly describe the aberrations. In this paper, 
the phase is expanded on a truncated Zernike basis. The impact 
of using a regularization metric with such a basis is studied later 
in this paper. In the MAP framework, the regularization metrics 
1Z(4> U ) and lZ(<fi d ) are deduced from the assumed a priori statis- 
tics of 4> u and 4> d . Assuming these aberrations are zero-mean, 
Gaussian, and neglecting a priori correlations between Zernike 
modes, we obtain, for an estimation performed on N Zernike 
modes: 



K{<t> x 



- 2 a x R a 



1^ 



07 



2 ^ < 

fc=l Xk 



(4) 



where a Xk is the assumed phase variance per Zernike mode, R ak 
the covariance matrix, and a x a.N element vector containing the 
estimated Zernike coefficients a Xk . Here x is either u (upstream) 
or d (downstream). 

The minimization of metric J(a, (3, <fi u , <p d ) of Eq.(3) is per- 
formed by means of a limited memory variable metric (BFGS) 
method (Press et al. 2007; Thiebaut 2002), which is a fast quasi- 
Newton type minimization method. It uses both gradients J^- 

and Flux a and offset /3 are analytically obtained using 

gradients ^ and ^ (implementation details, including gradi- 
ent expressions, can be found in Appendix A). 

Sauvage et al. (2012) established that a suitable diversity 
phase <p div for COFFEE was a mix of defocus and astigmatism: 
$ dm = af v Z 4 + af v Z 5 with af v = af v = 80nmRMS, 
introduced upstream of the coronagraph. We therefore use this 
diversity phase in the following. 

2.2. Coronagraphic imaging model 

The imaging model used by COFFEE in the criterion minimiza- 
tion (equation (3)) requires a coronagraphic PSF expression. In 
this paper, we use the analytical coronagraphic imaging model 
developed by Sauvage et al. (2010), whose formalism is devel- 
oped in this section, where r is the pupil plane position vector, r 
its modulus, and 7 the focal plane position vector. The entrance 
pupil function P«(r) is such that: 



P u (r) = n 



with n ( 



(5) 



= 1 for r < 



pupil entrance diameter, 
otherwise, and $ is an apodization function. In this paper, 
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we consider that the impact of amplitude aberrations is negligi- 
ble, which is a reasonable assumption for a ground-based, high- 
contrast imaging system such as SPHERE. Considering only 
static aberrations (no residual turbulent aberrations), the electric 
field i&A in the entrance pupil plane can be written as 



Simulation 



*a0) = P u {r)e 



(6) 



The field amplitude ^3(7) in plane B can be calculated, fol- 
lowing Sauvage et al. (2010), using the analytical coronagraphic 
imaging model (which is called "perfect coronagraph model" 
hereafter): 



*b( 7 ) = FT-\^ A (r)) - TyoFT" 1 (P„(r)), 



where i]o is the scalar that minimizes the outcoming energy from 
focal plane B, whose analytical value is given by 



V0 = Jf Jj**A(r)Pu{r)dr, 



where 



TV = 



JJ s P* u (r)P u (r)dr. 



(8) 



(9) 



It is worthy mentioning that 770 is the exact definition of the 
instantaneous Strehl ratio given by Born & Wolf (1989). One can 
notice that ryo = 1 when there is no aberration upstream of the 
coronagraph (0 u (r) = 0), so that \J>£ = in such a case. No 
aberration in the entrance pupil leads to no outcoming energy 
from plane B, and thus to a perfect extinction in the detector 
plane D. 

Propagating the wave from plane B (Eq. (7)) to plane D, we can 
write the electric field ^jj(-f) in the detector plane: 



*d(7) =FT- 1 jp d (r)e J W > "M+ < W T '»} 



(10) 



where Pd(f) is the Lyot stop pupil function: Pd(r) = 

Tl(^j^j P u (r), with Dd the Lyot stop pupil diameter (Dd < 

D u ). For the sake of simplicity, we omit the spatial variables r 
and 7 in the following. The coronagraphic PSF of the instru- 
ment, denoted by h c , is the square modulus of *S?d - 



h c (0„,0 d )=|FT- 1 (P d e 



j(<t>«+<t>d) 



?7 oFr- 1 (P de ^<0| 



(11) 



In this paper, this expression of the coronagraphic PSF is the 
one used by COFFEE for estimating <f> u and <fi d ; i.e., Eq. (11) 
is inserted into the imaging model (Eq. (1)) used in the criterion 
minimization described in Eq. (3). 

As described by Sauvage et al. (2010), this model, which an- 
alytically describes the impact of a coronagraph in an imaging 
system, considers that the coronagraph removes the projection 
of the incoming electric field on an Airy pattern, represented by 
the parameter ?/o (Eq. (8)). Since it does not assume small aber- 
rations, it can be used for any wave-front error upstream of the 
coronagraph. The quality of the fit of this analytical imaging 
model with the ARPM coronagraph is discussed later in this pa- 
per (Section 3.5). 



image size 

Light spectrum 

Aberration upstream of the 
coronagraph (4> u ) 
Aberration downstream of 
the coronagraph (4> d ) 
Zernike basis used for <f> u 
and 4> d simulation 



93x93 % (128x128 pixels, 
oversampling factor: 1.38) 
monochromatic (A = 635 
nm) 

WFE = 80 nm RMS 
WFE = 20 nm RMS 
36 Zernike polynomials 



Phase estimation: COFFEE 



Zernike basis used for <fi u 



(7) and <p d reconstruction 



36 Zernike polynomials 



Regularization metric 



none 



Table 1. COFFEE: simulation parameters used for the performance as- 
sessments of sections 3.1, 3.2 and 3.3 



3. Performance assessment by numerical 
simulation 

The aim of this section is to quantify the impact of each er- 
ror source on COFFEE'S aberration estimation. Such a study 
will show COFFEE's sensitivity to the classical error sources 
that limit the phase retrieval in a real system (and thus the fi- 
nal extinction of the coronagraph), which will be of high inter- 
est in defining COFFEE's upgrades. Likewise, it will allow us 
to estimate the accuracy level expected on our AO bench. In 
this section, we present the evolution of this reconstruction error 
with respect to the incoming flux (Section 3.1), to the size of the 
source (Section 3.2), to an error made on the assumed diversity 
phase used in the reconstruction (Section 3.3), and to the num- 
ber of Zernike modes used in the reconstruction (Section 3.4). 
For each error source, coronagraphic images will be computed 
using the imaging model presented in Eq. (1), using the per- 
fect coronagraph model to calculate the coronagraphic PSF h c 
whose expression is given Eq. (11). COFFEE will then perform 
the phase estimation using these two images. The compatibility 
of COFFEE with realistic coronagraphic images will be studied 
as well (Section 3.5) by computing coronagraphic images us- 
ing a realistic coronagraph model and then running COFFEE to 
estimate the aberrations both upstream and downstream of the 
coronagraph. 

Table 1 gathers the parameters used for these simulations. 

The chosen wave-front error (WFE) values upstream and 
downstream of the coronagraph for these simulations are typi- 
cal of the aberrations that will be estimated on our AO bench in 
Section 4 (so that experimental results can be compared to the 
following simulations). Since these simulations are performed 
with a small number of Zernike modes (36), there is no need of 
regularization metrics in such simulations. 
To simulate realistic aberrations, we have considered that the 
variance per Zernike mode er^ was decreasing with the radial 
order n(k) of the considered Zernike mode k (Noll 1976): 



1 



n(ky 



(12) 



This corresponds to a decrease in the static aberration spatial 
spectrum as j^p, where v is the spatial frequency, which is a 
common assumption for mirror fabrication errors. To evaluate 
COFFEE's performance, we define the reconstruction error e x 
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Flux (ph.) 
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Fig. 2. Aberrations upstream (<f> u (WFE = 80 nm), top) and down- 
stream (4> d (WFE = 20 nm), bottom) of the coronagraph: reconstruc- 
tion error (solid red line) as a function of the incoming flux a. For 
comparison, i (cyan dashed line) and (magenta dashed line) theo- 
retical behaviours are plotted for detector noise only and photon noise 
only (respectively). 

(x stands for u (upstream) or d (downstream)) as 



iV-l 

e = a ^ |ofc - a k \ 2 (13) 

\ fc=2 

with a,k the Zemike coefficients (starting with k = 2 correspond- 
ing to tilt) used for the simulation, Sft the reconstructed Zernike 
coefficients, and N the number of Zernike modes. In this sec- 
tion, every reconstruction error value is an average value, com- 
puted from ten independent simulated phases. 

3. 1. Noise propagation 

The ultimate limitation of an instrument lies in the amount of 
noise in the images. In Figure 2, we present the reconstruction 
error for the aberrations upstream (</>„) and downstream (4> d ) of 
the coronagraph with respect to the total incoming flux. Photon 
noise and detector noise (<7d et = 6 e _ ) are added in the corona- 
graphic images for simulation. 

The evolution of the reconstruction error presented in Fig- 
ure 2 is proportional to (1/a) for the detector noise limited 
regime (low flux) and to (1/y/a) for the photon noise limited 
regime (high flux). In this figure, it can be seen that for an in- 
coming flux a > 10 6 photons, the reconstruction error e u for the 



0.0 0.1 0.2 0.3 0.4 0.5 

FWHM (X/D) 

Fig. 3. Error reconstructions upstream (red line) and downstream 
(blue line) of the coronagraph as functions of the size of the source on 
the coronagraph. 

phase upstream of the coronagraph is smaller than 1 nm RMS. 
Thus, in a calibration process, where high values of flux (> 10 6 
photons) can be easily reached, COFFEE'S performance will not 
be significantly affected by noise. 

It is noteworthy that the results of many similar simulations with 
various levels of upstream aberrations show that COFFEE'S re- 
construction error does not depend on the amplitude of the aber- 
rations upstream of the coronagraph, as long as the diversity 
phase amplitude is larger than the WFE of the aberrations to be 
estimated. 

3.2. Impact of the source size on the reconstruction error 

Our imaging model, presented in Section 2.1 (Eq. 1), assumes 
an unresolved object. Thus, the presence of a real source with a 
given spatial extension will have an impact on the phase recon- 
struction, which is quantified here. We consider here a Gaussian- 
shaped laser source, emitted from a single-mode fiber. Because 
of the incoming light coherence, it can be represented as a Gaus- 
sian amplitude in the entrance pupil plane (where COFFEE as- 
sumes a uniform amplitude). Knowing this, coronagraphic im- 
ages are simulated by considering a small coherent Gaussian- 
shaped beam (FWHM < 0.5-^) on the coronagraph, and then 
processed by COFFEE. 

Since the imaging model assumes an unresolved object, both 
reconstruction errors for the phases upstream and downstre am of 
the coronagraph increase with the FWHM of the coherent object, 
as showed in Figure 3, but remains low: for an FWHM smaller 
than , the reconstruction error is indeed sub-nanometric. The 
size of the laser source will thus definitely not be a limitation 
for COFFEE: if this error is not negligible in the total error bud- 
get, it is possible to include it in the imaging model used by 
COFFEE (Eq. 1) as a non-uniform (Gaussian) entrance pupil 
function P u (r). 

3.3. Sensitivity to a diversity phase error 

The diversity phase <f> div = af^Z^ + af lv Z 5 has been de- 
fined in Section 2.1. This phase <p div is one of the inputs that 
COFFEE needs in order to perform phase retrieval, so it must 
be calibrated as accurately as possible. To optimize the use of 
COFFEE, the impact of an error on such a calibration is stud- 
ied. In this section, we consider that the diversity phase used to 
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Simulation 



Fig. 4. Error reconstructions upstream (solid red line) and downstream 
(solid blue line) of the coronagraph as functions of the error on the 
diversity phase. 



create the diversity image is not perfectly known. The corona- 
graphic simulated diversity image is computed with a diversity 
phase (f>' cUv = 4> dlv + 4> err , with 4> err a randomly generated 
phase of given RMS value, and COFFEE'S phase reconstruction 
is done considering that the diversity phase is equal to <f> div . In 
Figure 4, we see that the reconstruction error increases linearly 
with the calibration error on the diversity phase, with a slope of 
0.5. Thus, the requirement on the calibration precision for the di- 
versity phase is typically the precision wanted for the aberration 
measurement. 



3.4. Impact of aliasing 

The phase estimation is performed here on a truncated Zernike 
basis. In real images (recorded from a bench), some speckles 
will originate in high-order aberrations. These aberrations, 
which cannot be fitted by the truncated Zernike basis, will 
have an impact on the phase estimation, called aliasing error 
hereafter. Thus, it is necessary to study this aliasing error as 
a function of the number of Zernike modes used in the phase 
reconstruction. Here, we generate a phase on a large number of 
Zernike modes, and compute the corresponding images using 
the perfect coronagraph model. Aberrations both upstream 
and downstream of the coronagraph are then estimated by 
COFFEE using an increasing number of Zernike modes. Since 
one of the aims of this simulation is to determine the size 
of the truncated Zernike basis to be used with experimental 
data recorded on an in-house bench, the noise level in the 
simulated images corresponds to the one we have on this 
bench. The total incoming flux is 5 10 6 photons, and the 
detector noise is cr<j et = le" per pixel. Parameters used for 
this simulation are gathered in Table 2. This simulation has 
been done with and without a regularization metric, so that we 
can demonstrate the relevance of this metric on phase estimation. 

Figure 5 presents the evolution of the reconstruction errors 
when the number of reconstructed Zernike modes increases. 
Here, every reconstruction error (Eq. (13)) is calculated on a 
basis of 350 Zernike modes; thus, the error originates both in 
high-order aberrations, which are not considered by COFFEE 
because of the Zernike basis finite size (modelling error), and in 
the impact of these high-order aberrations on the estimated ones 
(aliasing). The WFE corresponding to the aberrations that are 



image size 

Light spectrum 

Aberration upstream of the 
coronagraph (4> u ) 
Aberration downstream of 
the coronagraph (4> d ) 
Zernike basis used for <p u 
and <f> d simulation 
Incoming flux 

noise 



93x93 % (128 x 128 pixels, 
oversampling factor: 1.38) 
monochromatic (A = 635 
nm) 

WFE = 80 nm RMS 

WFE = 20 nm RMS 

350 Zernike polynomials 
5 10 6 photons 

photon noise, detector noise 
(Odet = 1 e~) 



COFFEE: phase estimation 



Zernike basis used for <p u 
and <fr d reconstruction 
Regularization metric 



from 15 to 275 Zernike 

polynomials 

With and without 



Table 2. COFFEE: simulation parameters for studying the aliasing er- 
ror. 



not estimated by COFFEE (from N to 350, where N varies be- 
tween 15 and 275 according to Table 2) is called "unmodelled 
WFE" hereafter. 

In the plot of the reconstruction error upstream of the corona- 
graph evolution (Figure 5, top), one can see that without a reg- 
ularization metric, the reconstruction error increases for a large 
number of Zernike modes. An interpretation of this behaviour 
is the following: because high-order aberrations have a smaller 
variance, their associated speckle intensity is lower. Thus, ow- 
ing to the photon and detector noise in the image, the SNR is 
smaller for these aberrations. Such behaviour leads to a trade-off 
between aliasing and noise amplification for the optimal number 
of Zernike modes (Figure 5). The best number of Zernike modes 
is then a function of the aberrations level (WFE) and spectrum, 
as well as of the level of noise. The use of a regularization met- 
ric allows us to avoid this noise amplification (Figure 5): the re- 
construction error roughly reaches a saturation level (rather than 
growing to very high values). Additionally, the use of regular- 
ization reduces the aliasing error, and avoids the need for the dif- 
ficult and somewhat ad hoc choice of number of Zernike modes 
for the reconstruction. 

According to the results presented in Figure 5, we have chosen to 
estimate the aberrations upstream and downstream of the coron- 
agraph on 170 Zernike modes with the regularization metric of 
Eq. (4). 

3.5. Model mismatch 

We have already demonstrated that ARPM images are com- 
patible with the perfect coronagraph model and therefore with 
COFFEE estimation in Sauvage et al. (2012). The Roddier & 
Roddier Phase Mask (RRPM) (Roddier & Roddier 1997; Guy on 
et al. 1999) consists in a 7r phase shifting mask slightly smaller 
than the Airy disk. Additionally, the use of a circular prolate 
function as entrance pupil apodization <I>p (ARPM), proposed 
by Soummer et al. (2003), leads in a perfect case (no aberrations 
upstream of the coronagraph) to a total suppression of signal 
in the detector plane. In the simulations presented hereafter, 
realistic ARPM coronagraphic images are computed following 
Soummer et al. (2007) to consider an accurate numerical rep- 
resentation of Lyot-style coronagraphs. Then, we use COFFEE 
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50 100 150 200 250 

Number of reconstructed Zernike modes 




o L i i i i i i i i i i i i i i i i i i i i i i i i l_ 

50 100 150 200 250 

Number of reconstructed Zernike modes 

Fig. 5. Error reconstructions upstream (top) and downstream (bottom) 
of the coronagraph as functions of the number of reconstructed Zernike 
modes, with a regularization metric (solid blue line) and without (solid 
red line) 




20 40 60 80 

Upslreom WFE {nm RMS) 



Fig. 6. Error reconstruction upstream of the coronagraph with respect 
to the WFE of the aberration upstream of the coronagraph 



to reconstruct both phases upstream and downstream of the 
coronagraph. Here, when using the formalism developed in 
Section 2.2, the prolate apodization function <I>p is included in 
both simulation and reconstruction imaging models. 

Because the perfect coronagraph model is not exactly 
identical to an ARPM (although their responses to aberrations 




Fig. 7. Adaptive optics testbed schematic representation. M; : fold mir- 
rors; MP;: parabolic mirrors; L;: lenses (doublets); BS: beam splitter; 
TTM: Tip-Tilt mirror; DM: Deformable mirror; RRPM: coronagraphic 
focal plane mask; <I>: prolate apodizer; WFS: AO wave-front sensor 

is very close), there is a model mismatch in the estimation 
of aberrations upstream of the coronagraph <fi u , which varies 
linearly with the WFE of <p u , as shown in Figure 6. The model 
mismatch can thus be quantified as 7.5% of the WFE RMS 
value of 4> u , except for very small WFE (< 1 nm RMS), where 
the variation is non-linear, but remains below 1 nm RMS. 
Since the variation in this model mismatch varies linearly with 
the WFE of <p u , it should not limit the ability to compensate for 
the aberration upstream of an ARPM using COFFEE as focal 
plane wave-front sensor (WFS). 



4. Laboratory demonstration 

In this section we present experimental validations in the corona- 
graphic phase diversity. These validations are done on the bench 
BOA, described in Section 4. 1 . Section 4.2 describes a carefully 
designed method developed to introduce calibrated static aberra- 
tions on the AO bench to be measured with COFFEE. The error 
made on the measurements of aberrations upstream of the coron- 
agraph (NCPA) is quantified in Section 4.3. Section 4.4 presents 
the static aberration measurement performance, and Section 4.5 
details the procedure for compensating for the measured aberra- 
tions. 

4.1. Experimental setup 

Figure 7 shows the design of our in-house bench. The input 
beam, emitted from a fibered laser source (A = 635 nm) comes 
through the prolate apodizer which is in the entrance pupil 
plane (P u )- The beam is reflected by the tip-tilt mirror (TTM) 
and then by the deformable mirror (DM, entrance pupil, D u = 
40 mm, 6x6 actuators). The beam-splitter sends a fraction of 
the beam to the AO wave-front sensor (Shack-Hartmann, 5x5 
sub-apertures). On the other channel, the light is focused onto 
a RRPM, whose diameter is d c = 18.1 /im (angular diameter 
is 1.06-^y-). After going through the Lyot stop plane (Pd, with 
Dd = 0.99D U ), the beam is focused onto the camera (256 x 
256 pixels images with an oversampling of 2.75, detector noise 
Cdet = 1 e~). For faster computations, recorded images are re- 
binned to 128 x 128 pixels images with an oversampling of 1.38. 
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- DM 

3. The half difference <fi cal 



is our estimate of <j> 



cal " 



Fig. 8. Introduction of calibrated aberration on BOA: case of a pure 
spherical aberration. Left: theoretical wave-front (top) and DM intro- 
duced wave-front (bottom). Right: corresponding Zernike modes for 
the theoretical introduced aberration (solid red line) and the DM intro- 
duced aberration (dashed blue line). 

4.2. Introduction of calibrated aberrations 

To evaluate COFFEE'S performance, we introduce calibrated 
aberrations on the bench using a process described in this sec- 
tion. We consider an aberration phase <p cal to be introduced on 
BOA. First, since the phase is represented by the DM with a fi- 
nite number of actuators (6 x 6), the introduced aberration will 
not match the aberration <fi cal perfectly, as illustrated in Figure 8 
in the case of a pure spherical aberration. 

Our aim is here to introduce, using the DM, the closest 
aberration to the aberration 4>cai- We let F be the DM influ- 
ence matrix (obtained by calibration); any DM introduced aber- 
ration <p DAI can be described as a set of actuator voltages u 
(4> DM = Fu). We are thus looking for the set u ca i which solves 
the least-squares problem: 



U ca l 



argmin||FM- (p cal 



The solution of this problem can be written as 



Ucal 



Tcf> c 



(14) 



(15) 



with T the generalized inverse of matrix F. Using the inter- 
action matrix D (resulting from calibration), we can compute 
the corresponding set of slopes s ca i (s ca i = Du ca i), which can 
then be used to modify the AO loop reference slopes s le f . Thus, 
closing the AO loop with the reference slopes s le f + s ca i, we in- 
troduce an aberration 4>®Jf = Fu ca i = FT(f) cal on the bench, 
which is the best fit of <fi cal in the least squares sense. 
We also have to consider that the bench BOA presents its own 
unknown static aberrations 0^ OA and <^ 0A upstream and down- 
stream of the coronagraph (respectively). Thus, if a calibrated 
aberration <p cal is introduced in the entrance pupil, aberrations 
4> u upstream of the coronagraph will be 



<f>u = <t>cal + <t>l° A - 



(16) 



To get rid of the unknown aberration 0^ OA , we perform a differ- 
ential phase estimation: 

1. We introduce the aberration (j)^ 1 on the bench. A phase 

- + - DM ~ BOA 

4> u = 4> C ai + 4> u I s estimated using focused and diverse 
images recorded on the camera. 

2. The opposite aberration —(fi^ 1 is then introduced. A phase 

• DM - BOA 



4>u = -4>cai + 4>u is estimated. 



The first use of this process is to calibrate the diversity phase 
itself. Since this phase will be introduced using the AO system, 
the actually introduced diversity phase will not exactly match 
the theoretical mix of defocus and astigmatism. We introduce 
the aberrations <p div and —4>div on tne t> encn using the AO sys- 
tem. These two aberrations are then estimated using classical 
phase diversity (no coronagraph), with a pure defocus of diver- 
sity phase introduced using a flat glass plate of known thickness 
e in a focused beam. 

Such a process gives us an accurate estimation of the diversity 
phase really introduced on the bench, with an estimated accuracy 
of 4 nm RMS on the introduced aberration. This calibration is 
then used in COFFEE'S estimations performed on experimental 
images. 

4.3. Performance assessment: error budget 

From simulations presented in Section 3, we establish an error 
budget for estimating aberrations upstream of the coronagraph 
using experimental data: 

o Photon and detector noise error: on the BOA bench, the typ- 
ical incoming flux is /boa = 5 10 6 photons. Knowing that 
we have photon noise and a detector noise with = 1 e _ , 
we can evaluate the noise error: e no j se = 0.9 nm RMS. 

o The diversity phase 4> div has been calibrated using classical 
phase diversity, using the process presented in Section 4.2. 
Such an estimation has been performed with an error of 4.0 
nm RMS (value calculated from an error budget evaluated 
for a classical phase diversity estimation on the BOA bench. 
Such accuracy has already been obtained on this bench by 
Sauvage et al. (2007)). According to Section 3.3, this error 
on the diversity phase leads to an error e mo d e i = 2.0 nm RMS. 

o The source is a coherent Gaussian-shaped beam whose 
FWHM is 0.27^ on the coronagraph. According to the 
simulations of Section 3.2, this leads to a reconstruction 
error: e D bj = 0.7 nm RMS. 

o Residual turbulent speckles, which originate in uncorrected 
turbulent aberrations, are not included in the imaging 
model. To measure the impact of these speckle on the 
reconstruction, several wave-fronts have been successively 
recorded using a commercial Shack-Hartmann wave-front 
sensor. From these acquisitions, we calculate the WFE of 
the residual turbulent phase: = 1.2 nm RMS. This 

residual turbulence will create speckles on the detector, 
which will be considered by COFFEE as originating in 
NCPA. Thus, the residual turbulence error e tul -b made by 
COFFEE is estimated to e tur b = = 1.2 nm RMS. 

o Aliasing error, which originates in high-order aberrations, 
has been studied in Section 3.4. For a phase upstream of 
the coronagraph estimated on A*" = 170 Zernike modes, we 
have easing = 18.3 nm RMS. 

o From simulations, we know that the model mismatch is 7.5% 
of WFE. For this study, we will not estimate aberrations with 
a WFE stronger than 80 nm RMS. For such a WFE, the 
model error is e mo d e i = 6.0 nm RMS. 

As one can see in Table 3, the error budget is mainly driven 
by the aliasing error. The second most important term is the 
model mismatch (even though it goes to zero with the WFE). 
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Error budget 



Noise 


^noise 


= 0.9 nm RMS 


Model mismatch 


^model 


= 6.0 nm RMS 


Error on diversity 


£div = 


2.0 nm RMS 


Resolved object 


£obj = 


0.7 nm RMS 


Residual turbulence 


Eturb = 


= 1.2 nm RMS 


Aliasing 


^aliasin 


, = 18.3 nm RMS 


Total error 


e = 
RMS 


VEi3 = 20.6 nm 



Total error per Zernike e' = 1.6 nm RMS per esti- 
mode mated Zernike mode 



Table 3. COFFEE: error budget for the estimation of an aberration up- 
stream of the coronagraph on BOA. 



4.4. Measurement of aberrations upstream of the 
coronagraph 

In this section, we introduce calibrated aberrations on the BOA 
bench upstream of the coronagraph, and then estimate them with 
COFFEE in order to evaluate its performance. In the course of 
this study, we realized that the position of the coronagraphic im- 
age on the detector (quantified by the tip-tilt downstream of the 
coronagraph) is a critical issue. Indeed, it occurred that COFFEE 
was able to perform phase retrieval only for downstream tip- tilt 
[a 2 , a 3 ] values within the range [-100 nm RMS; 100 nm RMS] 
([— el?])- To 8 et rid of this constraint, we have developed a 
method to perform a preliminary estimation of the tip-tilt down- 
stream of the coronagraph. This method, which uses the diver- 
sity image, is fully described in Appendix B. 

4.4.1 . Measurement of tip-tilt upstream of the 
coronagraph 

We present the estimation of a tilt aberration upstream of the 
coronagraph using COFFEE in this section. Using the AO sys- 
tem, we introduce a tilt aberration by adding a constant value 
dsrr to the AO wave-front sensor references slopes s re f, and then 
closing the AO loop on the slopes s re f + Ssjj. To accurately 
calibrate the introduced tilt, for each position, we first estimate 
the aberrations using classical phase diversity (no coronagraph). 
Then, the RRPM is put in the focal plane, and the same opera- 
tion is repeated: for each position, we record two images, and 
then estimate the aberrations using COFFEE. 

From the upstream tilt reconstruction performed by COF- 
FEE (Figure 9), we calculate an average reconstruction error: 
e t ii t = 2.1 nm. Part of this error is due to an error on the estima- 
tion of tip-tilt downstream of the coronagraph. An improved es- 
timation has been performed by setting boundaries on the down- 
stream tip-tilt. Its value is evaluated before COFFEE'S estima- 
tion using the method described in appendix B with the diversity 
coronagraphic image recorded for a tip-tilt upstream the coron- 
agraph value close to nm RMS (centered coronagraph). Such 
an estimation process gives us an estimation of tip-tilt down- 
stream of the coronagraph {a^ , af'} with an accuracy of ±1.5 
nm RMS. Using this estimation as the starting value for the min- 
imization, and setting bounds of ±1.5 nm RMS on it, we pro- 
cessed the same experimental data. This, in turn, results in a 
better estimation of tilt upstream of the coronagraph (Figure 9), 
with an average error e t iit = 1.5 nm, which is close to the ex- 
pected error per Zernike mode given in Section 4.3 (e' = 1.6 nm 
RMS). 
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Fig. 9. Estimation of a tilt aberration on BOA: calibration (solid blue 
line) and COFFEE'S estimation with bound on the tip-tilt downstream 
of the coronagraph (dashed crossed red line) and without boundaries 
(dashed diamond green line) 

4.4.2. NCPA measurements 

In this section, we introduce aberrations upstream of the coron- 
agraph. The aberration <fi cal is expanded on the first 15 Zernike 
modes (which is the largest number of modes we can properly 
describe with our 6x6 DM), and then we estimate these aberra- 
tions using COFFEE, following the process described in Section 
4.2. To take the DM action into account on the introduced phase 
(illustrated in Figure 8), aberrations <p cal are first estimated with 
classical phase diversity (no phase mask in the coronagraphic 
focal plane (Sauvage et al. 2007)). This estimation gives us a 
calibration of the introduced aberration, which is then used to 
evaluate the accuracy of COFFEE'S estimation. 

At convergence of the reconstruction, a very good match can 
be observed between the experimental images and the ones com- 
puted for the estimated aberrations (Figure 10, top and middle). 
This, in turn, results in a very good match between the aberra- 
tions measured by COFFEE (Figure 10, right) and the introduced 
ones (Figure 10, left). 

From the experimental phase estimation presented in Figure 10, 
we compute a reconstruction error between the classical diver- 
sity phase calibrated aberration and COFFEE'S estimation: 

e exp = 22.5 nm RMS. (17) 

One can notice that this error is close to the expected error bud- 
get, i.e. that there is a good match between the performance 
assessment study carried out in Section 3 and the experimental 
results presented in this section. 

4.5. Low-order NCPA compensation 

Lastly, the ability of COFFEE to compensate for the aberrations 
upstream of the coronagraph is experimented on BOA. In Sec- 
tion 4.4, the aberrations upstream of the coronagraph are ex- 
panded on 170 Zernike modes, in order to have the smallest re- 
construction error (according to Section 3.4). 
As previously mentioned, the compensation on BOA is limited 
to the 15 th Zernike mode. Thus, what is required in a closed 
loop process is the most accurate estimation of 15 Zernike modes 
rather than an accurate measurement of every estimated Zernike 
mode. Using a basis of 36 Zernike modes for the reconstruc- 
tion is sufficient to give an accurate estimation of the first 15 
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Fig. 10. COFFEE: NCPA estimation of an introduced phase (j> cal 
on BOA. Top: for an aberration +4> ca i> recorded coronagraphic image 
from the bench (left) and computed image using the reconstructed aber- 
ration tp u (right) (log. scale, same range for both images). Middle: 
same images for an aberration —4> cal introduced and a reconstructed 

aberration <j> u (log. scale, same range for both images). Bottom: cal- 
ibrated introduced aberration (left) and COFFEE estimated aberration 
(right) 

Zernike modes: the aliasing error, which is the most important 
error source, will mainly degrade the estimation accuracy of the 
reconstructed high orders (close to Zse)- 

To demonstrate the ability of COFFEE to be used in a closed 
loop, we introduce a set of aberrations on the DM by modify- 
ing the reference slopes, as described in Section 4.2. Then, we 
use the pseudo-closed loop (PCL) method described in Sauvage 
et al. (2007). This iterative process has two stages: for the PCL 
iteration i: 

1. acquisition of the focused v c and diverse v, images; 

2. estimation of the aberration cj> u upstream of the coronagraph; 

3. computation of the corresponding reference slopes correc- 

tion 5s = gDT<p u , where D and T are the interaction and 
influence matrices defined in Section 4.2 and g is the PCL 
gain; 

4. the AO loop is closed on the modified reference slopes. 

The computation time (step 2) varies from 1 minute to 2.5 
minutes, allowing us to compensate for quasi-static aberrations 
upstream of the coronagraph. This compensation process is lim- 
ited by the estimation accuracy of the first 15 Zernike modes 
performed by COFFEE, which corresponds to the error budget 
established in Section 4.3), and by the ability of the DM to re- 
produce a given wave-front. Indeed, the correction introduced 



Pseudo-closed loop : Voridnce 




0-0 
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Fig. 11. PCL on the bench BOA (<?pcl = 0.5): variance of the resid- 
ual static aberrations upstream of the coronagraph for the 36 COFFEE 
estimated Zernike modes (solid red line) and the 15 corrected modes 
(solid blue line). The magenta dashed line represents the ultimate per- 
formance one can reach according to the error budget detailed in 4.3 

on the bench (step 2 of the PCL process) is the best fit of the 

a i 

estimated phase <f> u in the least-square sense (as presented in 
Section 4.2). The difference between the estimated aberration 
and the actual introduced correction will thus limit the compen- 
sation performance of the PCL process. Considering these two 
limitations, one can compute the variance ct| a (f° r tne ^ TSt 15 
Zernike modes) that can be reached on the BOA bench: 

cr| 0A = 4.4 10~ 2 rad RMS 2 . (18) 

The correction and stabilization of the NCPA variance can 
be seen in Figure 11. One can see that the variance of the 15 
corrected Zernike modes reaches the expected asymptotic value 
°boa- This result is the very first demonstration of COFFEE'S 
ability to compensate for aberrations upstream of the corona- 
graph. A compensation at levels compatible with SPHERE or 
GPI-like instruments will require th using a DM with many more 
actuators, and working on the reduction of the dominant term of 
the error budget, which is aliasing. 

5. Conclusion 

In this paper, we have presented a thorough simulation study 
(Section 3) and a first experimental validation (Section 4) of the 
coronagraphic wave-front sensor called COFFEE, which con- 
sists mainly in the extension of the phase diversity concept to a 
coronagraphic imaging system. From the validation and careful 
performance assessment of COFFEE, we showed that COFFEE 
is currently limited by the aliasing error, due to high-order aber- 
rations, which are difficult to model with a Zernike basis. 
In Section 4, we presented a first experimental validation of 
COFFEE using an ARPM. We introduced calibrated aberrations 
upstream of the coronagraph (NCPA), using the AO sub-system, 
and estimated them with COFFEE. The accuracy we obtained on 
these estimation shows a very good match with our error budget. 
Lastly, we used COFFEE in an iterative process to perform a 
preliminary validation of COFFEE'S ability to compensate for 
the aberrations upstream of the coronagraph. 
Several perspectives are currently considered to optimize COF- 
FEE: firstly, in order to minimize the impact of the aliasing error 
on the phase reconstruction, we plan to perform the phase re- 
construction on a pixel-wise map, which is more suitable than a 
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truncated Zernike basis. Secondly, we would like to improve the 
imaging model, both to make COFFEE work with other coro- 
nagraph than the ARPM and to reduce the model error, which 
is currently the second most important one, even though it goes 
to zero with the WFE. Two solutions are considered. In the ab- 
sence of residual turbulence, an accurate imaging model is ob- 
tained by propagating the electric field through each plane of the 
coronagraphic imaging system (Figure 1) for an arbitrary focal 
plane coronagraphic mask. Such a method, where no model er- 
ror needs to be considered, can be used for a laboratory calibra- 
tion. Alternatively, a more accurate analytical imaging model, 
which could include a residual turbulent aberration, can be de- 
veloped. Such a model, which could include a residual turbulent 
aberration, will ultimately allow us to perform NCPA estimation 
on images from the sky. These improvements should allow us 
to estimate and compensate for the aberrations upstream of the 
coronagraph using COFFEE with a nanometric precision in a 
closed loop process. 

A further perspective is to extend COFFEE to phase and ampli- 
tude aberration estimation, in order to create a dark hole region 
in the coronagraphic image. 
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Appendix A: Implementation details 

COFFEE performs a phase estimation by minimizing a criterion 
J whose expression is given by equation 3. To estimate <fi u and 
4> d (expanded on a truncated Zernike basis), we need both gra- 
dients and where a x = {a Xl , a X2 , a XN } is a vector 
that contains the Zernike coefficients, for an aberration expanded 
on N Zernike modes (x is for u (upstream) or d (downstream)). 
Let us write the numerical expression of J foc , using the notations 
defined in Section 2.1: 



t 



< OC [i]-«.h d et[<]*hH*]-/3 



foe r 



r foc 



C[t]-a.h dei [t}*h™[t]-P 



(A.l) 



= J 



foe 



-K 
J 



4> d 

div i 



K 



<t> d - 



With t the pixel position in the detector plane. er° c and er° IV are 
the noise variance maps. Considering the expression of J, we 
derive J foc , and then deduce the gradients expressions of J dlv us- 
ing a trivial substitution. Expressions of the regularization terms 

gradients are given by 



8a x 



(A.2) 



The calculation of gradients J^- and J^- is done following 

Mugnier et al. (2001): first, we calculate the gradient of ,P with 
respect to the PSF h c : 



83 



foe 



dtit 



[ah det (a.h dQt * h 



foe 



•foe \ 



(A.3) 



Then, the calculation consists in derivating the gradient of 
the PSF h c with respect to phases <p u [k] and <p d [I] at pixels k, I 
in pupils upstream and downstream of the coronagraph, respec- 
tively, and applying the chain rule, as already done in a non- 
coronagraphic case, e.g. in Thiebaut & Conan (1995). The cal- 

culation of both gradients and 



dJ" 



d<t> d [l] 



gives 
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dc/> u [k] 
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8 h f ° c V 



d<t> u [k\ 4- o h f r 



E 



t 



dh 



foe 
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foe 



FT 



djt { T 



[I] 



(A.4) 



(A.5) 



With 3 and 9f? the imaginary and real part (respectively), and 



dr]o 



jP 2 .e jq 



V(<K, <t> d ) = P d e^ + ^ *(cf> u ,4> d ) = FT" 1 ^) 
M<i> d ) = Ptf?* d * d (4> d ) = FT~ 1 (ip d ). 

(A.6) 

Since the phases are expanded on a Zernike basis, we need 
the gradients of J foc with respect to the Zernike coefficients a Xi 
of phase <p x . These gradients are given by the expression (Mug- 
nier et al. 2001): 

8 P oc 

(A.7) 



da Xi 



E#^ 



84> x 



Flux a and constant background j3 are also analytically esti- 
mated during the minimization, considering that 

2 



-iP[t] + a.h d et[i]*fcP[i] + /3 



r P n[t] 



(A.8) 



Where p is for "foe" (focused) or "div" (diverse). For the sake 
of simplicity, we shall omit the variable t. We have 

8JP 

da 



E(/ldet * h V c ) 2 ^ /l det * h P c 
P — + ^2^ — ^~ 



E 



„P 2 

er „ 



(A.9) 



CTn 

h det ★ hi 



dJP _ y> 
Which gives us, in a matricial form: 



1 i p 
^E 3~2 ~~ E 3~2 



(h dc ,*hP) 2 



h<tet*W 




(A. 10) 

A simple matrix inversion gives us the analytical estimation 
of the flux a and the background /3 for each iteration. 
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begin phase reconstruction. This method performs, on our ex- 
perimental images, a fast preliminary estimation (~ 1 second 
for a 256 x 256 image) of the tip-tilt downstream of the coron- 
agraph with an accuracy of 1.5 nm RMS, which is far enough, 
compared to the level of accuracy (±100 nm RMS) required by 
COFFEE. 



Fig. B.l. Coronagraphic diversity images computed for an aberration 
4> u + 4>div upstream, cf> d downstream of (left) and the only diversity 
aberration 4> div (right). The shape of both images is mainly driven by 
diversity aberration. 

Appendix B: Tip-tilt estimation downstream of 
the coronagraph 

The tip-tilt downstream of the coronagraph (which represents the 
image position on the detector) strongly limits COFFEE'S per- 
formance. Indeed, we determine that the phase estimation was 
accurate when -100 nm RMS < a, < 100 nm RMS, with <n the 
Zernike coefficient for tip or tilt (i <G {2, 3}). Beyond this range, 
COFFEE is unable to properly estimate both phases <p u and <p d . 
Such a phenomenon strongly limits COFFEE'S performance on 
a bench, since its utilization requires a restrictive location of the 
PSF on the detector. 

To get rid of this limitation, we have developed a simple and 
fast method of estimating the tip-tilt downstream of the coro- 
nagraph before COFFEE'S estimation, based on the diversity 
image. This image is created by adding a known aberration 
<t> div = a f v Z 4 + af v Z 5 (af v = af v = 80 nm RMS) to 4> u . 
Since the amplitude of this aberration is important (a<p div =113 
nm RMS), the speckles we have in the coronagraphic diversity 
image mainly originate in this diversity aberration. This is il- 
lustrated in Figure B.l, where we show two diversity images: 
one computed with randomly generated phases <p u (WFE 30 nm 
RMS), 4> d (WFE 10 nm RMS), and another computed with no 
aberrations other than the diversity ones. 

As one can see in Figure B.l, we can clearly identify the 
aberrations which originate in the diversity 4> div . The principle 
of our method lies in the research of these well-known aberra- 
tions (since we know the phase 4> div we introduce) in the diver- 
sity image i A by comparing it with a theoretical diversity image 
r c , calculated with no other aberrations than the diversity ones: 

i d Cth = htex*h c (4> div ,4> d = 0). (B.l) 

The comparison of i A Cth with % & c is performed using the 
method developed by Gratadour et al. (2005), which consists in 
minimizing the following criterion Jtt 
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Jtt{x, y) 



•div / 



(x , Vo) - i c " h ( x o, Vo) * S{x -x,y a - y) 



div 



(B.2) 



Where 8 is the dirac function. Minimization of Jtt gives us 
the shift [xmi yj\/]between both images. It is then possible to 
calculate the corresponding tip (0,-2) and tilt (03) downstream of 
the coronagraph knowing the image sampling s: 



a 2 



2s XM 



03 



(B.3) 



Finally, these estimated tip-tilt values are given to COFFEE as 
an input of the minimization, and are used as initial values to 
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